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Abstract. We analyze the one-dimensional (ID) and the two-dimensional (2D) repulsive 
Hubbard models (HM) for densities slightly away from half-filling through the behavior of two 
central quantities of a system: the uniform charge and spin susceptibilities. We point out that 
a consistent renormalization group treatment of them can only be achieved within a two-loop 
approach or beyond. In the ID HM, we show that this scheme reproduces correctly the metallic 
behavior given by the well-known Luttinger liquid fixed-point result. Then, we use the same 
approach to deal with the more complicated 2D HM. In this case, we are able to show that 
both uniform susceptibilities become suppressed for moderate interaction parameters as one 
take the system towards the Fermi surface. Therefore, this result adds further support to the 
interpretation that those systems are in fact insulating spin liquids. Later, we perform the same 
calculations in 2D using the conventional random phase approximation, and establish clearly a 
comparison between the two schemes. 



1. Introduction 

The anomalous electronic properties of the metallic phase of the copper-oxide high-Tc 
superconductors are widely believed to be linked to the manifestation of a non-Fermi liquid 
state arising from the two-dimensional (2D) character of those compounds. Motivated by this, 
physicists have then been trying to understand theoretically what are the precise conditions in 
which Landau's Fermi liquid paradigm ^[21 (i-e. the scheme based on the concept of low-lying 
electronic quasiparticles in the system) breaks down in 2D models. 

There are a number of approaches for accomplishing that task in 2D systems. The most 
direct route is through Landau's theory of spontaneous breaking of a continuous symmetry, 
often induced by an instability with respect to a relevant external perturbation. That case 
usually leads to a phase transition scenario accompanied by the realization of a long-range 
ordered state in the system. Moreover, it implies the onset of gapless low-energy elementary 
excitations: the so-called Goldstone modes (see Goldstone's theorem in Refs. |31 ID)- This 
theory provides the framework for understanding, among several other examples, the long-range 
magnetic order realized in the Mott insulating state of the 2D Hubbard model (HM) exactly at 
half- filling, and the resulting emergence of gapless antiferromagnetic magnon excitations, which 
manifest themselves solely through the spontaneous breaking of the spin rotation symmetry in 
the system. 



Alternatively, another viable approach for that is produced quantum-mechanically by the 
disordering induced by strong fluctuations resulting from low-dimensionality effects themselves 
or by frustration effects in 2D models. These interacting regimes in turn do not have any 
broken symmetry realized in the ground state and, as a result, all the correlations are purely 
short-ranged in real space. In view of this, in principle, there can be no gapless excitations in 
those cases, since Goldstone's theorem does not apply anymore. These exotic states cannot be 
described by the conventional Landau's theory mentioned before, and are commonly referred 
to as spin liquid ground states. Their existence was first predicted by Anderson long ago in 
the context of a Heisenberg model defined on a 2D triangular lattice jS], and in [Hj a decade 
later, and they might possibly hold the key for the appropriate theoretical explanation of the 
underlying mechanism of the high-Tc superconductors. 

To understand these issues better, theorists have been generally looking for similar physics 
in simpler models. For instance, in strictly one-dimensional (ID) systems, the situation is far 
more understood in that respect [71|Hl|nj. It is by now very well-established that Landau's Fermi 
liquid description breaks down completely in those particular cases. In fact, the appropriate 
effective low-energy theory for the ID metallic state is instead the so-called Luttinger liquid. 
This happens due to specific features associated with ID, i.e, the always nested Fermi surface 
(FS) consisting of only two discrete points {+kp and —kp), and the very restricted phase space 
available for the elementary excitations above the corresponding ground state. Consequently, 
the low-energy excitation spectrum becomes fundamentally dominated by bosonic collective 
excitations given by two types: the holons and the spinous excitations. These particles in turn 
imply an exact separation of the charge and spin degrees of freedom in the system. As a result, 
they are most effective in destroying completely the fundamental picture of well-defined low-lying 
quasiparticles, which are the basis of Fermi liquid theory. 

In this paper, it is our main goal to investigate the low-energy properties of the well-known 
2D repulsive HM for densities slightly away from half- filling. To do this, we analyze two central 
quantities of the model: the uniform charge and spin susceptibilities. We choose these two 
physical quantities in view of the fact that they naturally provide important information about 
the existence or not of low-lying excitations in the system. As will become clear soon, to perform 
a full renormalization group (RG) calculation of those quantities, it is essential to set up at least 
a two-loop order RG calculation for it. As a testing ground of this approach, we first calculate 
these quantities in the ID repulsive HM away from half-filling. In this case, we show that both 
quantities approach finite values, and the resulting low-energy description is indeed a metallic 




Figure 1. The Feynman diagrams for the one-particle irreducible four-point function 
r^^HPi)P2)P3)P4) in the backscattering channel up to two-loop order (third order). The free 
propagators are represented by either solid or dashed lines according to their association with 
the corresponding branches. 



phase (in fact a Luttinger liquid) in agreement with well-known exact results (Bethe ansatz, 
bosonization, etc). Later, we apply the same considerations to the more complicated 2D HM, 
which is the case we are mainly interested here. We point out that there is an interaction regime, 
where both quantities become strongly suppressed as one approaches the low-energy limit, thus 
adding further support to the interpretation that the resulting effective theory should be given 
by an insulating spin hquid state Then, we compare these results with the conventional 
random phase approximation (RPA) case, and show that this latter approximation is not able 
to detect clearly such an exotic ground state in the system. Finally, we discuss the physical 
origins of the discrepancies of those two results. 



2. Testing ground: The ID repulsive Hubbard model 

First of all, we define the Hubbard model (HM) defined on an ID lattice here. This model 
describes a system of tight-binding electrons hopping only to nearest neighbor sites, and 
interacting mutually through a local (in real space) repulsive interaction parametrized by C/ > 0. 
The corresponding Hamiltonian given in momentum space is 



where the single-particle energy dispersion is simply = —2t cos ka — /i, and ip].^ and V'fco- are 
the usual creation and annihilation operators of electrons with momentum k and spin projection 
a =1, |. Besides, fi stands for the chemical potential, a is the lattice spacing, and NsUes 
is the total number of lattice sites. Another important parameter here is the width of the 
noninteracting band, which is given hy W = 4:t. 

The FS of the noninteracting system is given by two discrete points i+kp and —kp). Since 
we are mostly interested in the low-energy properties of this model, we linearize the energy 
dispersion about the two Fermi points, i.e., = (|^| — kp), where vp = {d(,k/dk)\k=kp is 
the Fermi velocity of the system. We restrict the allowed single-particle states to lie within the 
interval kp — X < \k\ < kp + A, where A is a fixed ultraviolet (UV) microscopic momentum 
cut-off, and is the corresponding energy cutoff given hy = 2vpX. As a result, we are 
able to identify two different regions in k-space: the branch linearized about k = +kp, 
and the "-" branch linearized about k = —kp. Thus, the corresponding free propagators may 
be represented by two different lines in the Feynman diagrams: solid line for the former, and 
dashed line for the latter. In addition, regarding the interaction term, we use the standard g- 
ology parametrization j7j , namely, the interaction processes involving large momentum transfers 
(near 2kp) are described by the gi coupling, whereas the processes involving small momentum 
transfers between particles located at different branches are described by the 52 coupling. The 
former interactions are usually called backscattering processes, while the latter are referred to 
as forward scattering processes. In what follows, we will deliberately neglect the contributions 
due to the so-called Umklapp processes, which are described by the coupling (73, in view of 
the fact that we will not consider here the exceptional case of an exact half-filled band in 
the model. Moreover, we will make one further approximation and neglect forward scattering 
pocesses between particles belonging to the same branch, which are conventionally described by 
the coupling constant (74. 

In the following, we will use the so-called field-theoretical RG strategy to tackle this problem. 
In view of this, it is more convenient to work here with the Lagrangian of the model instead of 
the Hamiltonian. Thus, we write down the Lagrangian of the ID HM as 
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Figure 2. The self-energy diagrams up to two-loop order. The one-loop order (first order) 
diagrams are the nonsingular Hartree-Fock terms. The first nonanalyticity comes from the 
two-loop order (second order) diagrams. 



k,a,a=± 

13 [92Sa56/3f - giSa-f^lSs] {P + Q - k,t) V'{_)^ {k, t) V'(-)/3 {q, t) V'(+)a {p, t) , 

(2) 

where ipj^^^ and are now fermionic fields associated to electrons located at the ± branches. 
The summation over momenta must be appropriately understood as J2p = Wl^^r) / dp in the 
thermodynamic limit, where V is the "volume" of the ID system. Since this should represent 
the ID HM, the Lagrangian is written in a manifestly SU{2) invariant form and, in addition to 
this, the couplings are initially defined as gi = §2 = iy / sites)U ■ 

Since the HM is microscopic, the couplings and the fermionic fields in the Lagrangian are 
defined at a scale of a few lattice spacings in real space (i.e at the UV cutoff scale in momentum 
space). These parameters are usually inaccessible to every day experiments, since the latter 
probe only the low-energy dynamics of the system. In RG theory, these unobserved quantities 
are known as the bare parameters. In fact, when one tries to construct naive perturbative 
calculations with such quantities, one is immediately confronted with several Fcynman diagrams, 
which turn out to be singular in the low-energy limit. These are the so-called infrared (IR) 
divergences in field theory. They appear in both backscattering and forward scattering four- 
point vertex corrections (see, for instance. Fig. 1), and in the calculation of the self-energy 
alike (Fig. 2). Those divergences do not imply necessarily that perturbation theory is simply 
unable to give any reliable prediction of the low-energy properties of the model, but only that 
the way it is formulated is not the appropriate one for this case. Therefore, one must generalize 
the perturbative approach so as to eliminate such singularities in the calculations. This is the 
strategy of the so-called renormalized perturbation theory, i.e, one rewrites all unobserved bare 
parameters in terms of the associated renormalized (or observable) quantities. The difference 
between them will be given by the so-called counterterms, whose essential role is to absorb 
all divergences Tip to a given order in perturbation theory. If this program is successfully 
accomplished, then the theory is said to be properly renormalized. 

Following this, we have 



gi ^Z-"^ {giR + AgiR) , (z = 1 , 2) , 
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Figure 3. The Feynman diagrams for the uniform response function T^'^'^\p,q ~ 0) up to 
two-loop order. The triangular vertices stand for the renormalized response function, whereas 
the diagram with a cross represents the counterterm. 



where '^(a)a- ^^'^ respectively, the renormalized fields and couplings, whereas Aip^-^^ and 

AgiFt are the corresponding counterterms. Besides, Z is the quasiparticle weight and measures 
the coherence of the quasiparticle picture in the effective (renormalized) description. As we 
explained in our previous paper jllj . this parameter is related to the self-energy effects, and as 
such is essential for a consistent two-loop renormalization group analysis of the model. 

Nevertheless, there are still many ways to define the appropriate counterterms using the above 
prescription. In order to solve this ambiguity, we must establish a renormalization condition for 
the observable parameters. We use here the one-particle irreducible four-point and two-point 
functions (for standard definitions see, for instance, Ref. ^2]). Therefore, we have 

r^''^(p01 +P02 =P03 -POl = =P3 = kF,P2 = PA = -kp) = -%ij(w), (5) 

r^'^^(poi +P02 =P03 -Poi = uj,pi = Pa = kF,P2 = Ps = -kp) = -ig2R{uj), (6) 

and 

r^^\ko = uj,k = kF)=io, (7) 

where the subscript "0" just means that we are dealing with the energy component associated 
with the external legs. Besides, uj is an energy scale which denotes the proximity of the 
renormalized theory to the FS. In this way, to flow towards the FS, we let a; —> 0. 

To compute the RG flow of the renormalized couplings, one needs to consider that the original 
theory (i.e. the bare theory) does not know anything about the extrinsic scale lo. In other words, 
the bare couplings are independent of uj. Thus, using the RG condition udgi/duj = 0, we obtain 

^ dgiR ^ gjji ^ gfn 

duo TTVp 2tt'^Vp 

^dg2R ^ gin ^ gfn 
duj 2'KVp Air'^Vp 



(8) 
(9) 



In these equations, we note that there are two possibihties of flow as one approaches the FS. 
First, for models defined with an initially repulsive backscattering coupling (i.e. gi > 0), the 
low-energy effective theory naturally scales to a line of fixed points given by g^j^ = and 
92R 7^ 0- Therefore, these systems arc in the same universality class of the so-called Tomonaga- 
Luttinger model where backscattering processes are neglected altogether. For this reason, they 
are often referred to as Luttinger liquids. This is the case for the ID HM, in which we are 
mainly interested here since it will provide a testing ground for our later analysis of the 2D 
HM. We note, in passing, that there is another possibility of flow which happens when the 
initial model includes a backscattering coupling which is attractive (i.e. gi < 0). In this case, 
the backscattering coupling becomes relevant in the RG sense, and the resulting theory is in 
the same universality class of another well-known model in ID, which goes under the name of 
Luther-Emery model. For simplicity, we will only deal with the first case scenario in this paper. 

We next derive a flow equation to analyze the quasiparticle weight Z of the Luttinger liquid 
universality class as one approach the FS. Calculating the Feynman diagrams of Fig. 2, and 
using the renormalization condition of Eq. (7), we obtain 

= ^ (9m + 91r - 91R92R) = 7 (10) 

where 7 is the so-called anomalous dimension. In the vicinity of the fixed point g\j^ = and 
92R 7^ 0, the quasiparticle weight scales as Z ~ (a;/0)'>'*. Since 7* > 0, Z becomes suppressed at 
the FS if we let a; — . This result asserts that there are no coherent low-lying quasiparticles 

in these systems. 

In order to obtain the uniform susceptibilities of the system, we must first calculate the 
linear response function due to an infinitesimal uniform external field, which couples with the 
occupation number operator. Thus, we add to the Lagrangian the new term 

-hexternal ^«V'(a)a(^'' (1^) 

p,a=± 

where, again, everything is written in terms of the bare parameters of the model. This 
will generate an additional vertex (the one-particle irreducible uniform response function 
r(2,i)(p^^ ~ 0)), which will in turn contain new IR divergent diagrams but only at two- loop 
order or beyond in perturbation theory. As we will point out later in the paper, at one-loop 
order, there are no singular diagrams present in this calculation. Therefore, these contributions 
do not need to be regularized and, as a result, do not enter in the RG flow equations. For this 
reason, they are usually dealt with by other means such as the RPA approximation. On the 
other hand, at two-loop order, one can verify the existence of several divergent diagrams in the 
calculation of the uniform response functions (see Fig. 3). This important result will, in fact, 
allow a direct computation of their associated RG flow equations, which will eventually describe 
how these parameters vary as we take the physical system towards the FS of the system. 

Analogously, we rewrite the bare uniform response function in terms of its renormalized 
counterpart and the appropriate counterterm AT^ in the following way 



(12) 



As we pointed out above, one still needs to establish how the experimentally observable response 
are to be defined. We choose the canonical renormalization condition, i.e., T^'^'^^ijpQ = u,p = 
kp, g ~ 0) = —iT^^uj), where the scale lj plays the same role explained earlier. 

Now, we can define the two different types of uniform response functions, which simply arise 
from a symmetrization of the object with respect to the spin projection a, namely 



^cWM = T^^H+T^^(a;), (13) 
TilM = T^^(a;)-T^^(a;), (14) 

where T^f^^^g^ and are the uniform response functions in the charge and spin channels, 

respectively. To derive the corresponding flow equations, we recall that udTa/duj = 0. Thus 

.d^R _ ( 9ir92r \ , . 

.d^R _ ( 9ir92r\ ^r , . 

We observe here that the RG flow equations associated with the charge and spin uniform 
susceptibilities are identical in form. This has a simple interpretation, i.e, there is no separation 
of the charge and spin degrees of freedom in our current model. In other words, the holons and 
the spinous elementary excitations known to exist in these systems have the same velocity of 
propagation here. This result can be traced back to our initial Lagrangian, where we left out 
from our analysis the so-called 5(4 processes associated with forward scattering processes between 
particles from the same branch. These interaction processes are interesting in their own right, but 
they are widely recognized to be extremely difficult to be incorporated consistently in a full RG 
scheme [7]. Nevertheless, as we will see next, it is still possible to extract important information 
from our model concerning the low-energy properties of the Luttinger liquid universality class. 

As we saw earlier, the Luttinger liquid fixed point is defined by the Tomonaga-Luttinger 
model, i.e., 5*^ = and ^2/? 7^ 0- Due to the limitations associated with a two- loop order 
truncation in the present RG scheme, such a fixed point is reached only logarithmically in the 
RG fiow. On the other hand, it can be proven by other means, i.e. by using the so-called Ward 
identities associated with an exact conservation of left and right particles, that this fixed point 
indeed holds to all orders of perturbation theory [13j . For this reason, this result is fundamentally 
nonperturbative in character. Thus, in order to properly take the physical system towards the 
Luttinger liquid regime in our approach, we must set from the start the backscattering coupling 
giR equal to zero in Eqs. (|T^ and (fTT)|) . Therefore, we obtain: 

i^-^r^harge = (17) 



'^^.'^sLn = 0, (18) 



d_ 

This means that both uniform response functions satisfy trivially the RG invariance condition at 
the Luttinger liquid fixed point. In addition to that, the corresponding uniform susceptibilities 
(see also Fig. 4) are given by 

1 r 1 2 

Xcharge(Spin)i^) = 2tTVf '•'^-^ " ^^^^ 

Thus, since the uniform response functions are RG invariant quantities at the Luttinger liquid 
fixed point, so will be the corresponding uniform susceptibilities of the system in the same 
regime. This is related to the fact that the low-energy effective theory indeed contains both 
gapless charge and spin collective excitations, which are indicative of a metallic ground state. 
In this way, we can conclude that the two-loop RG approach is able to reproduce, quite simply, 
some very important aspects of the so-called Luttinger liquids, namely, that its fixed point 
corresponds to a metallic state with no fermionic quasiparticles excitations, and the resulting 
dynamics being dominated exclusively by bosonic charge and spin collective excitations jTl |H1 El • 
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Figure 4. The Feynman diagram associated with the renormalized uniform charge and spin 
susceptibilities. 

3. The 2D repulsive Hubbard model 
3.1. The two-loop RG approach 

Our two-loop RG calculation worked quite well in reproducing some low-energy properties of 
the ID HM such as the gapless nature of the resulting in both charge and spin excitation 
spectra. This result encourages us to try to apply the same approach for problems, in which 
there are no known exact solutions to this present date. This is what happens very generally 
in strongly correlated 2D models, due to the fact that well-known analytical techniques such 
as the bosonization and the Bethe ansatz are not easily extendable to such higher-dimensional 
systems. Here, we discuss the application of the two-loop RG method developed earlier to the 
2D repulsive HM slightly away from half-filling. 

In 2D, the repulsive {U > 0) Hubbard Hamiltonian defined on a square lattice is given in 
momentum space by 



where the energy dispersion is given by = — 2i [cos(fca;a) -|- cos{kya)] — fi, and Vko- ^■iid tpka 
are the creation and annihilation operators of electrons with momentum \c = {kx,ky) and spin 
projection a =|, j. The other definitions follow those of the ID model, i.e., n stands for the 
chemical potential, NgUes is the number of total sites, and a is the spacing of the 2D lattice. A 
very important difference here is the width of the noninteracting band, which is now given by 



The electron band filling of the system is controlled by the ratio fi/t. When fi/t = the 
model is exactly at half-filling condition. As we start doping it with holes, n/t takes slightly 
negative values. As we can see from Fig. 5 (a), the FS of such lightly doped system is nearly flat, 
and contains no van Hove singularities in it. That FS will be our starting point for the two-loop 
RG approach. However, this 2D problem has additional complications compared to the ID case, 
since we need now two momenta kx and ky to prescribe a renormalization procedure towards the 
FS. To simplify the calculations, we can perform a simple rotation of the axes by 45" degrees, 
and change these two components by a component parallel and another perpendicular {kj_) 
to the FS. As a result of this, the 2D problem will bear a more transparent resemblance to the 
ID case (see Fig. 5 (b)), and we can use the machinery developed in the previous section with 
some minor adaptations. To emphasize further the similarities of both problems, we divide the 
2D FS into four different regions (two sets of solid and dashed line patches). Here, we restrict 
the momenta at the FS to the flat parts only. Analogously to the ID problem, the interaction 
processes connecting parallel patches of the FS are always logarithmically IR divergent due to 
quantum fluctuations. In contrast, those connecting perpendicular patches always remain finite, 
and do not contribute to the RG flow equations in our approach. For convenience, we restrict 
ourselves to one-electron states labeled by the momenta k\\ = kx and k± = ky associated with 
one of the two sets of perpendicular patches. The momenta parallel to the FS are restricted to 




(20) 



W = 8t. 




Figure 5. (a) The half-filled Fermi surface (FS) of the 2D Hubbard model (dashed line), the 
lightly hole-doped FS (solid line), and (b) the latter FS after a rotation of the axes by 45° 
degrees. 

the interval —A < fcy < A, with 2A being essentially the size of the flat patches. The energy 
dispersion of the single-particle states is given by = vp i\k±\ — kp), and depends only on the 
momenta perpendicular to the FS. We define a label a = it which refers to the flat sectors at 
k± = =b/ci?, respectively. In addition, we take fci? — A < \k±\ < kp + X, where A is the fixed 
ultraviolet (UV) microscopic momentum cut-off with Q = 2vf\. 

In the same spirit as before, we now write down the Lagrangian of the 2D HM as 

^= '0(a)a('^'*)[^^t-'^^(l^-L|-^i^)]V'(a)<T(k,t) 
k,(T,a=± 

~^ [52505(5/37 - 51 '5a7MV'(+)5 (P + q - k,t)V'(_)^ (k,t)V'(-)/3 (q, *)''/'(+)« (P>*) > 

p,q,k":/5.7.'5 

(21) 

where all the definitions follow the ID HM case with and V'(±) being fermionic fields 

associated to electrons located at the ± patches and the summation over momenta being 
appropriately understood as J2p = W(27i")^ / d?P in the thermodynamic limit. Since this should 
represent the 2D HM, the Lagrangian is also written in a manifestly SU{2) invariant form and 
the couplings are initially defined as gi = §2 = (Y / sites)U . In addition, we again neglect 
Umklapp processes here. This choice is related to the fact that we are mainly interested in 
dealing with the HM slightly away from half-filling and, therefore, our FS does not intersect the 
so-called Umklapp surface at any point (see again Fig. 5 (a)). 

Apart from its physical importance especially regarding the possible explanation of high-Tc 
superconductivity, this 2D model is also very interesting from an exclusive RG point of view. 
In fact, when one deals with this FS model, the counterterms needed to renormalize the theory 
turn out to be continuous functions of the momenta parallel to the FS, rather than being simply 
infinite shifts. As a result, the couplings initially defined as constants in the 2D Hubbard 
Lagrangian acquire a momentum dependence and become renormalized functions of the three 
parallel momenta when one takes the system towards the FS. Thus, the new prescriptions become 



^ia)a - + AVgy (fell) = (22) 

4 

5j ^ ]J Z"^/^(Pi||) 9ii?(Pl|hl'2||,P3||) + A5'i/?(Pl|hP2|hP3||) • (23) 
i=l 

where = piy +^211 ~ P3\\- before, there are still several ways to define the appropriate 
counterterms using the above prescription. In order to solve this, we establish renormalization 
conditions for the one-particle irreducible four-point functions as follows 

^^^\P1\\,P2\\,P3\\,P01+P02 =P03-P01 = ^,Pl =P3 = -P2 = -P4 = kp) = -igiR{Pl\\,P2\\,P3f^), 

(24) 

^^^\P1\\,P2\\,P3\\,P01+P02 =P03-P01 =^,Pl =P4 = -P2 = -P3 = kp) = -ig2R{Pl\\, P2\\, P3\\; ^) , 

(25) 

and also for the one-particle irreducible two-point function as 

T^^\ko = uj,k = kF,k\\) =u;, (26) 

where u is again the energy scale which denotes the proximity of the renormalized theory to the 
FS. 

Following the same procedure explained before, the corresponding RG flow equations for 
the couplings and quasiparticle weight of this 2D problem become coupled integro-differential 
equations, which are impossible to solve analytically. Therefore, one must solve them using a 
numerical approach. We choose here the standard fourth-order Runge-Kutta method. We refer 
the reader to our paper 11 for more details. 

When one considers moderate initial interactions, the corresponding quasiparticle weight Z 
becomes suppressed as one approaches the FS. This is a strong indicative that the resulting 
low-energy effective theory does not contain low-lying electronic quasiparticles in the system. 
In addition to this, instead of the renormalized coupling functions diverging as happens in well- 
known one-loop RG flows, they flow, for some particular choices of momenta, to strong coupling 
plateau values. This result means that those systems belong to the same universality class of a 
strongly coupled theory, whose detailed information must be determined from a consistent RG 
calculation of other quantities such as the uniform susceptibilities as we have done explicitly for 
the ID HM case. In fact, we will see here that the suppression of the quasiparticle weight of the 
system will produce interesting effects in the RG flows of those quantities '14'. 

First, we must calculate the linear response due to an infinitesimal external field perturbation 
as before. Again, we add to the Lagrangian the following term 

-hexternal '^Mi^^J)^ {p) tp(a)a (P) ^ (27) 

p,a=± 

This will generate an additional vertex (the one-particle irreducible uniform response function 
p(2,i)^p^ q 0)), which will in turn be afflicted by divergences given by the Feynman diagrams of 
Fig. 3 as in the ID case. Therefore, we rewrite the bare quantity 7^ in terms of its renormalized 
parameter 7^ and an appropriate counterterm A7^ as follows 



T,(P||) = Z-i(pij) [T„^(p||) + ATi?(pi,)J . (28) 
where the Z function is given by 



Z{p\\) = 1 - ^^^4^2 / dkiidqii [2g2R o g2R + 2giR o - gm o g2R - g2R o giR] In . (29) 

where the sign o means that we are omitting the momentum dependence of the renormahzed 
couphng functions. As was already mentioned, A7^ must absorb the divergences generated by 
the Feynman diagrams of Fig. 3. However, we may stih define that counterterm in several 
ways. To solve this ambiguity, we make a prescription establishing precisely that the is 
the experimentally observable response, i.e., 

,Po = ^jP± = kF;q ~ 0) = —iT^{p\\,uj), 

where oj is the RG energy scale parameter. 

Now we turn to the definition of the two different types of uniform response functions, which 
arise from the symmetrization with respect to the spin projection a, i.e. 

'^ChargeiPW,^) = r^^(p|| , + T^^(p|| , w), (30) 
'Ts'iiniPW,^) = r.^(p||,^)-T^^(p|i,^), (31) 



where T^^h^rge '^Spin ^^'^ uniform response functions associated with the charge and spin 
channels, respectively. To compute the RG flow equations of these response functions, one needs 
to recall that the bare parameters are independent of the renormalization scale u. Thus, using 
the RG condition ivdTa/duj = 0, we obtain 



_d_ 



TchargeiPw) = 32^4^2 ( / dk\\dq\\[2g2R o g2R + 2giR o giR + giR o g2R + g2R o giR 

+9lR O 92R + g2R ° giR - 2g2R O g2R - 2giR o - 2g2R o g2R - 2giR o giR\T^harge{l\\) 

+ j dk\\dq\\[2g2R o g2R + 2giR o g^ - gm o g2R - g2R o giR]T(?f^argeiP\\)^ > (32) 



-252/?, o g2R - 2giR o giji - 2g2R o g2R - 2giR o giR]Tsli^{q\\) + Ts^i„(p||) 

X J dk\\dq\\ [2g2R o g2R + 2giR o g^ - gm o g2R - g2R o giR]^ , (33) 

where the sign o represents again the momentum dependence of the renormalized coupling 
functions which were omitted for convenience. To find the solutions for these RG equations, we 
resort to the same numerical method as described before. 

Finally, once the response functions are obtained, we can calculate the flow of the uniform 
charge and spin susceptibilities of the system. They are given by the same diagram as in the 
ID case (Fig. 4). Therefore, we have 

R ^ f R 1^ 

Xcharge{Spin)i^) = 4^2^^ J.A [^Charge{Spin)iPh ^)\ ■ (34) 

Once more, we use the same numerical procedure to estimate those quantities. Our results are 
displayed in Fig. 6. We note in this plot that both uniform susceptibilities flow at the same 
rate as we approach the FS. As a result, again, there is no spin-charge separation effects in this 
model. One way to possibly circumvent this would be to add 54-type interaction contributions 
and work hard to implement an appropriate RG scheme in this system (for this question, see also 
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Figure 6. The RG flow of the uniform charge and spin susceptibiUties as we increase the bare 
interaction strength given hy U = {g/'irvF)t. Our RG step I is given by I = ln(fl/uj). 



Ref. for a discussion of a closely related problem in ID). Notice that, for moderate initial 
couplings with Z going to zero, the susceptibilities become strongly suppressed in the same low- 
energy limit in markedly contrast to the ID case, where they remain finite. Therefore, our result 
adds further support to the interpretation that the low-energy effective theory for those cases 
should be a fully gapped state in both charge and spin excitation spectra ^1]. Since such an 
effective theory has only gapful excitations present, it cannot be related to any spontaneously 
broken symmetry state. Consequently, it should stay in a liquid phase-like state associated 
with a short-range order in the system. These effects become even more transparent when we 
calculate the susceptibilities associated with several order parameters (see Ref. |16j). There, we 
show that all susceptibilities are finite as we vary the RG energy scale. Thus, it is not possible 
to any given order to become long-range correlated in this 2D system. This unusual quantum 
state predicted long ago by Anderson is commonly referred to as an insulating spin liquid (ISL), 
and it is widely believed to play a central role in the appropriate explanation of the underlying 
mechanism of high-Tc superconductivity. 



3.2. The RPA approach 

We now compare the two-loop RG result with another approach for calculating the uniform 
susceptibilities that appears widely in the literature ^1 ^j, i.e., the random phase 
approximation (RPA). This method is used for estimating those quantities at an one- loop RG 
level. At this order, the renormalized couplings are known to diverge at finite energy scales, below 
which the RG is simply unable to access. The approximation consists of the assumption that 
the remaining information is given by the one-loop diagrams of the uniform response functions 
(see Fig. 7). However, since these diagrams are not IR divergent, it is not possible to derive 
RG fiow equations for those quantities at this order. The resulting equations will be integral 
equations which must be solved self-consistently in strong contrast to our previous two-loop RG 
approach, which allowed us to derive proper flow equations. 

Let us consider Eq. (26). If one goes only up to one-loop order, the quasiparticle weight 
Z must be equal to unity, since one is not taking into account self-energy corrections in the 




Figure 7. The Feynman diagrams for the calculation of the uniform response functions within 
the RPA approach. 



calculations. In this way, we have 



r„(p||)=Ti^(p||) + AT«(p||). 



(35) 



Now, calculating the Feynman diagrams of Fig. 7 and establishing the same renormalization 
condition used in the previous section, i.e., T^"^'^^ {p\\,pQ = u;,p± = A;i?;q 0) = —iT^{p^^,Lu), we 
obtain the following expression 



dk\\ 



'^aik\\)9iRiP\\,k\\,P\\) - 13^0^(^11 )5'2i?(P||,A:||,fc|| 



(36) 



We note here that this equation is indeed not a RG flow equation. Thus, symmetrizing this 
quantity with respect to the spin projection a, i.e. 



'^ChargeiPW^^) 



rf{p\\,u;)+T^%\\,u;), (37) 
^5pm(P|h^) = r^''(P||,a;)-r^^(p||,a;), (38) 

we obtain for the uniform charge and spin response functions the following equations 



1 

'^ChargeiPf^) = '^Charge(jP\\) + 1^;^ j _^dk\\T^j,^^g^{k\\) 

X 9iR{p\\,k\\,P\\\^) - '^92R{p\\,k\\,k\\;uj) 



1 

Tspin{p\\) + J^^dkiiT^pi^ikii)giRipii,ki^,pii;Lu). 



(39) 



(40) 



which now must be calculated self-consistently. Here, we discretize the FS in exactly the same 
way as explained in the previous section. 

Once these quantities are estimated, we can calculate the uniform susceptibilities in the RPA 
approximation, i.e. 
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Figure 8. The RG flow of the uniform charge susceptibility as we increase the bare interaction 
strength given by C/ = {g/'KVF)t within the RPA approach. Our RG step I is given by 
I = ln(0/a;). 
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Figure 9. The RG flow of the uniform spin susceptibihty as we increase the bare interaction 
strength given by C/ = {g/'KVF)t within the RPA approach. Our RG step / is given by 
/ = \n{Vl/uj). 



X^ChargeiSpin) (^) = dp\\ [T^harge{Spin) iP\\ , ^)] ' • (41) 

The results for these quantities are displayed in Figs. 8 and 9. In the charge susceptibility 
plot, although the initial flow gives the impression that it becomes somewhat suppressed as we 
approach the FS, we observe in fact a strong upturn of the flow towards infinity in qualitative 
agreement with other numerical results |171 I18j . As for the uniform spin susceptibility, we 
note a tendency to flow to zero in the low-energy limit as we increase the initial coupling U. 
Therefore, even though one might argue that the RPA approach captures an initial tendency of 
both quantities to become suppressed, this is in fact rather approximate. As a result, one is not 
able to detect clearly the formation of the ISL state in this 2D model within this approximation. 
In contrast, our two- loop RG calculation provided a more solid evidence that the natural 



candidate for the low-energy effective description of the 2D HM shghtly away from half-fihing 
is indeed the ISL. The underlying reason for the different behavior of the two approaches lies in 
the following physical principle. Quantum fluctuation effects arc well-known to be very effective 
in destroying long-range correlations especially in low-dimcnsional systems. In fact, the more 
fluctuations are taken into account, the more likely long-range ordered states are transformed 
into short-range ordered ones. This also becomes clear as a result of our work. In accordance 
with this, we can infer that if we were to consider higher-order fluctuations such as a three-loop 
RG calculation or beyond, we would approach the ISL state even more rapidly in this 2D HM 
system. 



4. Conclusions 

Here we investigated the renormalization of the uniform charge and spin susceptibilities of the 2D 
HM slightly away from half-filling within a two-loop RG approach. In our calculations, wc took 
into account simultaneously both the renormalization of the couplings and self-energy effects. 
As a warm-up example, we first applied this approach to the ID HM away from half-filling 
and reproduced several features associated with the Luttinger liquid fixed-point: the irrelevance 
in the RG sense of the so-called backscattcring interaction processes, the complete absence of 
low- lying quasiparticles in the vicinity of the FS, and finally the finiteness of the both charge 
and spin uniform susceptibilities indicating that the resulting state is indeed metallic. 

Once the RG technique was explained in detail, we moved on to the problem we were 
mostly interested in, i.e., the 2D HM for densities slightly away from half-filling condition. We 
emphasized the strong similarities to the ID case, and adapted the technique to deal properly 
with this 2D model. We pointed out that, for moderate bare interaction U, the coupling functions 
flow to strong-coupling plateau values for particular choices of the external momenta. Besides, 
this regime is also characterized by a strongly suppressed quasiparticle weight Z, which is a 
strong indicative of the absence of low-lying electronic quasiparticles in the resulting low-energy 
effective theory. Then, in order to extract further information about this state, we calculated 
the uniform susceptibilities associated with the charge and spin degrees of freedom, and showed 
that both quantities renormalize to zero in the low-energy limit. This result implies that both 
charge and spin excitation spectra are fully gapped, and the resulting state should be given by 
an insulating spin liquid. 

In order to compare our result with other numerical estimates encountered in the literature, 
we derived a RPA approach to calculate the uniform susceptibilities within a one-loop RG 
scheme. As a result, we reproduced the well-known one- loop results, and showed that they were 
in fact not sufficient to obtain in a clear way the ISL low-energy effective theory in this case. 
This should be contrasted with the two-loop RG approach, where we do observe this state. We 
explained that the physical reason for the discrepancies of both results is that, in this particular 
2D problem, one must add more quantum fluctuations in order to effectively remove every single 
instability that could potentially develop into long-range order in the system. Therefore, if one 
could go beyond two-loop RG order, one would not change this picture qualitatively, and might 
well approach the ISL state even more rapidly. 
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